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1.  Introduction 


Low-frequency  airborne  radar  systems  are  a  promising  technology  for  the 
surveillance  of  foliage-obscured  stationary  and  moving  ground  targets  such  as 
humans  and  vehicles.'^  Although  much  effort  has  been  put  forth  in  this  area  over 
the  years,  many  unsolved  problems  remain  as  the  foliage  ecosystem  itself — with 
scatterers  of  various  shapes  and  sizes  in  multifarious  arrangements — is  inherently 
highly  cluttered,  and  therefore,  poses  an  intricate  challenge  for  target  detection  and 
recognition,  especially  for  a  sensing  platform  in  motion.  Determination  of  target 
detectability  ultimately  necessitates  a  thorough  fundamental  understanding  of  the 
target  and  foliage  scattering  phenomenology  (as  well  as  the  tree  clutter  statistics 
and  attenuation)  as  a  function  of  frequency,  polarization,  signal  incidence  angle, 
and  foliage  properties.  As  such,  a  full-wave  electromagnetic  solver  is  proposed  in 
this  study  for  scattering  and  propagation  characterization  in  the  foliage 
environment.  Specifically,  the  analysis  framework  employs  a  parallelized  finite- 
difference  time-domain  (FDTD)  algorithm  for  deriving  the  far-field  responses  of 
complex  realistic  tree  structures  at  P-band  (200-500  MHz):  the  backscattering 
behavior  of  a  single  tree  is  first  considered,  and  then  large-scale  simulations  are 
performed  to  examine  the  responses  of  a  forest  stand  in  both  the  frequency  and 
imaging  domains.  The  solver  is  intended  as  a  computational  tool  for  evaluating  the 
effectiveness  of  detection  modes  based  on  methods  such  as  airborne  synthetic 
aperture  radar  and  ground  moving  target  indication;  nevertheless,  the  results  and 
observations  featured  are  also  expected  to  be  useful  for  remote  sensing  studies  in 
which  the  primary  interest  lies  in  retrieving  forest  parameters  from  the  scattering 
return. 

2.  Simulation  Framework 


As  the  first  step,  3 -dimensional  representations  of  realistic  tree  structures  need  to 
be  systematically  constructed.  The  geometrical  models  in  this  work  are  created 
using  the  open-source  random  tree  generator  Arbaro,^  which  implements  the  rule- 
based  growth  algorithm  developed  by  Weber  and  Penn.^  Essentially,  each  tree 
structure  is  generated  by  following  a  series  of  recursive  rules  derived  from 
geometrical  observables.  The  tree  structure  in  its  entirety  (trunk,  multi-level  branch 
complex,  and  leaves)  can  be  exported  from  Arbaro  in  mesh  form.  Given  the  long 
wavelength  at  P-band,  the  effects  of  the  leaves  are  often  assumed  to  be  negligible; 
therefore,  only  the  trunk  and  branches  are  retained  here  for  the  electromagnetic 
simulations.  Shown  in  Figs.  1-3  are  some  tree  model  examples. 
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Fig.  1  Quaking  aspen  (Populus  tremuloides)  tree  model  with  increasing  structural  fidelity 
from  left  to  right 


Fig.  3  Eastern  cottonwood  {Populus  deltoides)  tree  model 

After  the  mesh  has  been  properly  processed  (through  operations  such  as  closing  of 
holes,  translation  and  scaling,  etc.),  voxelization  of  the  structure  is  performed  with 
the  parity  count  and  ray  stabbing  methods. The  scattering  responses  of  the 
resulting  voxel  grid  are  then  calculated  with  a  customized  far-field  FDTD 
algorithm.  The  implemented  code  is  fiilly  parallelized  and  runs  on  large 
distributed  computer  systems  using  the  message-passing  interface  (MPI) 
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framework.  In  effect,  the  simulation  domain  is  decomposed  into  rectangular  sub- 
domains  and  the  FDTD  equations  are  solved  separately  for  each  sub-domain  within 
one  MPI  process.  The  computational  approach  is  highly  scalable,  even  when  the 
number  of  MPI  processes  is  in  the  hundreds.  More  details  on  the  electromagnetic 
solver  can  be  found  in  our  papers. 

3.  Numerical  Experiments 

The  backscattering  return  of  a  standalone  tree  above  a  half  space  is  first  analyzed. 
The  tree  species  considered  is  the  quaking  aspen  (Populus  tremuloides),  as  shown 
in  Fig.  1 .  Both  the  tree  and  the  ground  are  assumed  to  be  electrically  homogeneous, 
with  relative  dielectric  constant  and  conductivity  ( ,  cr^ )  of  (13.9,  39  mS/m)  and 

(5.45,  20  mS/m),  respectively.  The  tree  has  a  height  of  7.4  m,  with  a  base  diameter 
of  34  cm  and  a  canopy  diameter  of  3.4  m.  The  co-polarized  responses  (radar  cross 
section  and  phase  difference)  derived  from  simulating  the  structure  with  different 
levels  of  geometrical  detail  are  displayed  in  Figs.  4-6,  as  a  function  of  frequency 
and  elevation  incidence  angle  6. .  The  results  are  obtained  by  averaging  the 

responses  over  36  azimuth  incidence  angles  <l)^ .  For  the  set  of  parameters  employed 

here,  it  is  seen  that  the  backscattering  mechanism  of  the  trunk-ground  dihedral 
provides  the  dominant  signal  response.  The  branch  canopy  has  a  tendency  to 
enhance  the  radar  cross  section;  this  is  more  evident  and  consistent  at  horizontal- 
horizontal  Qih)  than  at  vertical-vertical  (vv);  the  inclusion  of  the  secondary  branches 
in  this  case,  however,  does  not  significantly  affect  the  overall  return.  Also  note  that 
the  impact  of  the  branch  canopy  seems  to  decrease  as  the  incidence  angle 
approaches  grazing.  In  general,  the  vv  responses  on  average  are  lower — and  more 
sensitive  to  the  incidence  angle — ^than  the  hh  responses;  both  of  these  observations 
can  be  attributed  to  the  Brewster  angle  effect. 
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(a) 


(b) 


Fig.  4  Average  response  (co-polarized  radar  cross  section  and  phase  difference)  of  a  single 
tree  as  a  function  of  frequency  and  structural  fidelity  at  9-  =  45° .  Dotted  line  =  trunk;  dashed 
line  =  trunk  +  primary  branches;  solid  line  =  trunk  +  primary  and  secondary  branches. 
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Fig.  5  Average  response  (co-polarized  radar  cross  section  and  phase  difference)  of  a  single 
tree  as  a  function  of  frequency  and  structural  fidelity  at  0-  =  60°  .  Dotted  line  =  trunk;  dashed 
line  =  trunk  +  primary  branches;  solid  line  =  trunk  +  primary  and  secondary  branches. 
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Fig.  6  Average  response  (co-polarized  radar  cross  section  and  phase  difference)  of  a  single 
tree  as  a  function  of  frequency  and  structural  fidelity  at  =  75° .  Dotted  line  =  trunk;  dashed 
line  =  trunk  +  primary  branches;  solid  line  =  trunk  +  primary  and  secondary  branches. 

Next,  the  simulation  of  a  large  forest  scene  consisting  of  36  randomly  generated 
trees  is  considered  (Fig.  7).  All  the  trees  (with  their  full  branch  structures)  are 
assumed  to  be  of  the  same  species  as  above,  arranged  within  an  approximately 
circular  area.  The  dielectric  properties  are  as  previously  described.  The  FDTD 
computational  domain  has  dimensions  29  m  x  29  m  x  9  m  and  is  discretized  into 
1.9  billion  cells.  Parallelized  simulations  were  performed  on  a  Cray  XC30  system 
after  partitioning  the  scene  into  256  sub-domains.  Each  simulation  run  at  1  set  of 
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incidence  angles  required  413  CPU  hours.  The  average  response  of  the  forest  stand 
at  =  75°  is  shown  in  Fig.  8.  Imaging  of  the  scene  at  each  elevation  incidence 

angle  is  demonstrated  by  re- focusing  the  backscattered  fields  captured  over  a  50° 
angular  aperture  (from  =  -25°  to  25°,  in  2.5°  intervals).  As  seen  in 

Fig.  9  {6^=15° ),  the  locations  of  some  of  the  tree  trunks  can  be  inferred  from  the 

co-polarized  images.  The  images,  however,  are  also  populated  with  other  clutter 
effects  such  as  multipath  interactions  (i.e.,  coupling  among  trees)  and  speckle.  The 
results  at  =  45°  and  =  60°  are  similar  to  those  in  Fig.  9. 


Fig.  7  Forest  stand  with  36  randomly  generated  trees.  Average  dimensions:  tree  height  = 
7.3  m;  base  diameter  =  39  cm;  canopy  diameter  =  3.2  m. 
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(a) 


(b) 

Fig.  8  Average  response  of  scene  in  Fig.  7 
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(a) 


(b) 

Fig.  9  Co-polarized  images  of  scene  in  Fig.  7:  a)  vv;  and  b)  hh 
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4.  Conclusions 


Full-wave  simulations  of  scattering  from  trees  are  considered  in  this  study. 
Geometrically  realistic  representations  of  tree  structures  are  first  created  using  an 
open-source  random  tree  generator.  Mesh  processing  and  voxelization  of  the 
structures  are  then  carried  out  to  generate  grid  models  that  are  amendable  to 
numerical  analysis.  Finally,  far-field  electromagnetic  scattering  characterization  is 
performed  using  a  parallelized  FDTD  solver.  To  demonstrate  the  practicality  of  the 
analysis  framework,  the  responses  of  a  quaking  aspen  tree  are  investigated  as  a 
function  of  frequency,  polarization,  signal  incidence  angle,  and  structural  fidelity. 
Subsequently,  large-scale  simulations  are  undertaken  to  examine  the  responses  of 
a  forest  stand  in  both  the  frequency  and  imaging  domains. 
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